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Abstract: 

We suggest that the topological susceptibility in gluodynamics can be found in terms 
of the gluon condensate using renormalizability and heavy fermion representation of the 
anomaly. Analogous relations can be also obtained for other zero momentum correlation 
functions involving the topological density operator. Using these relations, we find the 9 
dependence of the condensates (GG), (GG) and of the partition function for small 6 and 
an arbitrary number of colors. 



1 Introduction 



The importance of exact results in QCD or pure Yang-Mills (YM) theory is hard to over- 
estimate. There is as yet no thorough understanding of their nonperturbative aspects 
which are of crucial importance in most instances where strong interactions are involved. 
In these circumstances, any exact model-independent statements (theorems) about non- 
perturbative properties of QCD or YM theory become especially valuable. It turns out 
that some exact results can be obtained even in the absence of detailed knowledge of the 
confinement and chiral symmetry breaking mechanisms. These theorems will survive any 
further development of the theory. On the other hand, they serve as a test case for any 
nonperturbative models and confinement scenarios, which should respect them. 

Inspired by a recent impressive breakthrough [|l| in supersymmetric (SUSY) theories 
(see e.g. P] for a review), one may wish to address the issue of available exact results 
for ordinary, non-supersymmetric QCD or YM theory which is not considered as a softly 
broken SUSY model, but is rather taken on its own. Only a few of them are known so far. 
One well known example is provided by the so-called t'Hooft consistency condition |^ 
which ensures the anomaly matching in terms of composite and building block particles. 
Another series of theorems is based on strict inequalities which yield a number of quali- 
tative, but rigorous statements for QCD 0]. Two more classes of exact results are given 
by anomalous Ward identities ^ and low energy theorems based on a combination of the 
Ward identities technique with renormalizability arguments 0. The power of the latter 
method is that it allows for exact evaluation of some nonperturbative observables in QCD, 
see e.g. for a review. Among the variety of results, obtained within this approach, 
exact calculations of zero momentum correlation functions of scalar gluon currents (to be 
discussed below) are of special interest for our purposes. 

In this paper we address zero momentum correlation functions of the topological den- 
sity operator a^/ (STr)^^ with itself and the scalar gluon current a^/ {87i)GG, in pure YM 
theory (gluodynamics)Q. We will make an attempt to find these correlation functions in 
terms of the gluon condensates {GG), (GG) without resorting to any particular approxi- 
mation, but rather appealing to the general properties of the theory such as renormaliz- 
ability and dimensional transmutation. In particular, we will argue that the topological 
susceptibility in gluodynamics can be found in terms of the gluon condensate within a 
particular regularization scheme. Although our arguments will be somewhat heuristic, 
there is a hope that our final Eq.(|2T|) is correct. This hope is supported by the fact that 
the same relation can be also obtained along with a very different line of reasoning |^ 
which implies the same regularization scheme. 

Our interest in the correlation functions of the above type is mostly related to the 
problem of the 6 dependence , which has been long discussed in the literature practically 
since the discovery of instantons [jlO[- We note that little is known about the 6 dependence 
in YM theory. Prior to a recent renewed interest in these issues in supersymmetric theories 
||TT| , |l], 1^, essentially all what was known about the 6 dependence in YM theory was the 
fact that the physics should depend on 9 through the combination 9/Nc. This conclusion 
follows from the analysis of the famous Witten-Veneziano mass relation for the t]' [|T^ 

^ We would like to recall that in QCD all zero momentum correlation functions of GG are suppressed 
in the chiral limit by the light quark masses. 
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obtained within the large approach. After the works of Witten and Veneziano, the 
discussion of the 9 dependence in the hterature has mostly switched to field theoretical 
models other than gluodynamics. One of the purposes of this letter is to find the 6 
dependence of lowest dimensional condensates and the partition function in YM theory 
for arbitrary number of colors Nc and any value of ^ -C vr. 

There exist many different reasons to study the correlation function of the above type 
and the 6 dependence of the lowest dimensional condensates. Here we would like to 
mention a few of them. First, such relations are needed for the construction of a low 



energy effective Lagrangian |jT3|] which can be used to study the vacuum structure of the 
theory. Such an effective Lagrangian allows one to find the 6 dependence for all (and not 
only for small) values of 6. Secondly, the latter results for the 9 dependence in YM theory 
and QCD can be used for the axion physics ^M. Finally, these results may be of interest 



for the lattice studies |]T5[ and nonperturbative models of the YM vacuum including, in 
particular, instanton liquid models ||T^. 

The starting point of our method is the well known relation 

(y ~ — 

^Gl^G^^^ = -hmjMmj,^)^, (1) 

where \E' stands for a fermion field of mass M. The meaning of this formula is that 
the topological density operator is equivalent to the fermion bilinear (0) in the external 
gluon field in the infinite mass limitQ. This observation determines our whole approach 
to the problem. For our purposes, we consider gluodynamics as the low energy limit of 
YM theory coupled to a heavy fermion fieldQ. The fermion mass is assumed to be much 
larger than the inverse confinement radius, which ensures that this heavy fermion does 
not affect infrared properties of the theory. As its mass can be made arbitrarily large, all 
calculations are under theoretical control, and no problem can arise with this definition. 
On the other hand, the inclusion of a heavy fermion requires a regularization on a much 
higher ultraviolet scale. For the latter, we will assume the Pauli-Villars regulators. We 
will therefore work in the double limit Aym <C M -C Mr, where M is the mass of the 
auxiliary fermion, M/j is the Pauli-Villars fermion mass, and Aym stands for the dimen- 
sional transmutation parameter in gluodynamics. We will argue that a combination of 
the Ward identities technique with the fermion representation of the anomaly, Eq.(0), and 
its analog for the scalar case together with renormalizability arguments allow one to find 
zero-momentum correlation functions of the topological density in gluodynamics in terms 
of the nonperturbative gluon condensates (GG) , {GG) , within a particular regularization 
scheme. 

Our presentation is organized as follows. In Sect. 2 we propose a method which relates 
the topological susceptibility in gluodynamics with the gluon condensate. In Sect. 3 we 
combine this result with the Witten- Veneziano |T^ approach to the U(l) problem at large 
Nc- The 7]' meson mass in the limit —>■ , Nc —>■ oo will be expressed in terms of the 
gluon condensate of pure YM theory. Sect. 4 deals with the problem of the 9 dependence 
in YM theory for small values 9 <^'k. Sect. 5 contains our conclusions. 

^Note that a different relation of the topological density with light flavour fermion bilinears has been 
used starting from in the study of the ^-dependence in QCD by current algebra methods. 
■^In what follows, we use the terms gluodynamics and YM theory interchangeably. 
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2 Correlation functions of topological density 



In this section we study the zero momentum two-point function of the topological density 
(topological susceptibility) in pure YM theory. In what follows we will use the short 
notation for zero momentum correlation functions 

{A,B) = limi / c/xe*P^'(0|T{A(x)5(0)}|0) (2) 

p-*0 J 

of arbitrary operators A[x) and B{x). In this notation, we want to evaluate the two-point 
function : 

P = (■^^ GnijGnu , - G auG ny) ■ (3) 

As is well known, this correlation function vanishes to all order of perturbation theory 
since the topological density operator is a total derivative. Therefore, it seems that 
explicitly nonperturbative methods (for example, the instanton approximation) are needed 
to evaluate the two-point function (|^). What will be argued below is that this correlation 
function can be analysed using general arguments appealing to renormalizability of YM 
theory, together with the fermion representation of the anomaly, Eq.(|I]). 

To proceed with our arguments, we would like to remind the reader how a similar 
problem of the scalar two-point function 

S = {z^G^yGfj_y , z^G^yG^y) (4) 

was addressed by Novikov et. al. (NSVZ) These authors have derived a fundamental 
low energy theorem fixing zero momentum correlation functions of scalar gluon currents : 

I j dx{Q\T{0{x) ^G\m^)yM = -g{0) (5) 

Here 0{x) is an arbitrary color singlet local operator of canonical dimension d and h = 
ll/3Nc stands for the first coefficient of the /3-function in pure Yang-Mills theory^. As a 
derivation of NSVZ theorem (|]) is rather simple, for the sake of completeness we would 
like to recall it here. One starts with a redefinition of the gluon field 

Gfj,u = QoGfj^u , (6) 

where Qq is the bare coupling constant defined at the cut-off scale Mr. Then the path 
integral representation immediately yields the relation 

I j dx{Q\T{0{x) G\Q)m = -^(Y;^(0) • (7) 

On the other hand, the renormalizability and dimensional transmutation phenomenon in 
a massless theory (either QCD with massless quarks or gluodynamics) ensure that 



(O) = const 



M^exp 



1 d 



The same formula holds in QCD in the chiral limit with the substitution h = ll/3A^c ~ 2/3ri/ where 
rij is a number of light flavors. 
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with the choice b = 11 /SN^ — 2 /Suf or b = 11 /SN^, respectively. Finally, performing the 
differentiation yields the final result (^. Using the NSVZ theorem (|^) for the diagonal 
case O = {as/12iT)G'^, which is proportional to the conformal anomaly — 6as/(87r)G^, we 
obtain 

One should note that relation (J^) implies the Wick type of the T-product in the correlation 
function @. Moreover, Eq.(|^) means that perturbative contributions are subtracted in 
both sides of Eq.(H). In this sense, the correlation function can be obtained by the 
differentiation of \og{Z/ZpT) in respect to I/qq (here Z and ZpT stand for the full and 
perturbatively defined partition functions, respectively, see Eq.(PT|) below )Q. 

We would now like to interpret Eq.(^) in a different way which makes use of the view 
of gluodynamics as a low energy limit of YM theory coupled to a heavy fermion of mass 
M 0. An idea similar to what follows was first implemented by Kiihn and Zakharov (KZ) 
to relate the proton matrix element of the topological density {N\GG\N) to the known 
quantity {N\G'^\N) fixed (in the chiral limit) by the conformal anomaly. A technical trick, 
which will be used below, was suggested in [|19| where, in particular, we reproduced by 
this method the KZ relation for the matrix element (A^IGGIA^). In the present work, a 
similar idea is applied to correlation functions of the topological density. 

Let us introduce chiral projections 

*L = ^(1 - 75)^ , = ^(1 + 75)^ (10) 

of the fermion field. We can now use the well known expansion of the scalar density 

- M (^i^^ + ^R^A ^ ^G^M,, . (11) 

Using (|Tl|), we put Eq.(||) in the form 

, a. 



,^^^G' , -M {^L^R + ^r^l)) = ^{-M {^L^R + ^r^l)) (12) 

The following observation is crucial for what follows. Generally, Eq.(|T2|) is just equivalent 
to d^), and one might therefore think that we would not obtain anything new in the 
course of this identical re- writing. However, one can note that, as long as the heavy 
fermion decouples, a "rotated" fermion should also do, because the physical content of 
the effective low energy theory can not be changed by a redefinition of a heavy fermion 
field (see the end of this section for more detail). In particular, one can consider the chiral 
rotations (change of variables) 

^l'^r exp{ia)^L-^R , ^R<ilL exp{-ia)^R-^L ■ (13) 



^ The same result was obtained in using canonical methods with Pauli-Villars regularization. 
To one loop order in regulator fields, it was found that perturbative contributions add the identity to 
both sides of Eq.(|). 

^ One should note that both the /?- function and the dimensional parameter Ay m become different when 
YM theory is coupled to a heavy fermion. However, the decoupling ensures that the relevant parameters 
of a corresponding low energy effective theory (i.e. pure gluodynamics) are its own /3-function and Aym- 
The problem of matching theories with different color or/and flavor groups has been discussed more than 
once in the literature, mostly in the content of supersymmetric models jl^, |20| . 
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Different transformation properties of tlie two terms in (^) under rotations (|T^ imply 
that Eq.(pr2D is equivalent to two relations 



a, ^2 



G' , M^L^n) = ^{M^l^r) , (14) 



■127r ' " 36 
, a 



, ^ G\M^r^l) = ^{M^r'^l) . (15) 

Now we can take the difference of these equations to extract the combination of interest. 
One can see that in this case the chiral fermion bilinears will combine into the pseudoscalar 
gluon density GO. The relevant expansion reads (cf. Eq.(|l])) 

- tM (^l'^r - ^r^l) ^ ^G,,G,, . (16) 
^ ^ on 

Using ([T6|), we transform (p!^Jl5| ) into 



i^G' , ^G,^G,^) - . (17) 



Normally, both sides of this equation vanish because of the wrong CP-parity. However, 
we may think of a "regularized" theory with a 5-term for ^ -C 1, where Eq.(|T^) would be 
perfectly sensible. In any case, such a "regularization" by the ^-terms is only needed at 
the intermediate step (0) of our derivation. The final answer will be valid for any value 
of 9, including 9 = 0. As a double-check of our procedure, one can note that Eq.(p!7D 
coincides with the general NSVZ formula (^) for the particular case O = {as/STT)GG. 

Now we proceed similarly with Eg . ([T7|) . We substitute the G^ term in the l.h.s. and 
GG term in the r.h.s. of ([T7|) by the heavy quark expansions ([Tl|) and (piBf), respectively. 
Then (|l^ becomes 

{^Gf^uG,, , M {^L^ilR + ^r^l)) = ^^{M {^l'^r - ^R^i)) (18) 
Using the identity 



and transformation properties under chiral rotations (|1^) , we see that ( p^Sf ) is again equiv- 
alent to two equations 

a ~ - 8 - 

{^G,,G,, , M^l^r) = ^-(Mf'i75^fi) , (19) 

OTT OO 

n ~ - 8 - 

{-^G,,G,, , M^^^i) = ^-(M^^75^l) . (20) 
87r 3b 

It is obvious that the scalar combination of Eqs. (|T^,^) exactly reproduces Eq. (|T^ which, 
as we have seen, is a particular version of the original NSVZ theorem (^. New informa- 
tion is contained in the pseudoscalar part of Eqs.(|l^,^). Taking the difference of these 
equations and multiplying the whole expression by i, we finally arrive at the relation 

,/,.(0|r{^G<5.^Ge}|0>^-lAG^>. (21) 
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The relation obtained is clearly scheme dependent. In fact, we have used above the 
particular scheme in which all condensates and correlation functions are defined through 
the path integral. Thus Eq.(^T]) implies the Wick type of the T-product, i.e. the two-point 
function (pT]) is the second derivative of log Z in respect to 0. Therefore, Eq.(pTD does not 
contain ultraviolet divergences which are present in the factor Zpx (see Eg. (pTD below) 
and drop out after the differentiation in 9. This ensures that the nonperturbative gluon 
condensate in Eqs. @ and (pTf ) is the same quantity. For more comments on the scheme 
dependence of our results, see the end of Sect. 4. We note that after Eq.(^) is established, 
arbitrary n-point correlation functions of operators G^, GG can be found using Eqs. (P), 
(p!7|) and (|2TD by further differentiation in respect to l/g^ and 6. 

Furthermore, as the whole line of reasoning leading to Eq.(^) may appear rather 
vague, we would now like to make further comments justifying all steps made in the 
derivation. First, let us discuss the use of the dimensional transmutation arguments (^) 
in our case where we introduce an auxiliary heavy fermion field. It might seem that 
in this situation Eq.(H) does not hold, since now we have the mass parameter M in the 
theory. Let us, however, notice that for a heavy fermion mass M any dependence on M in 
formulas like (|), besides a perturbative redefinition of Ay^^, can only enter via corrections 
in powers of Ayj^f/M which die off asymptotically as M cxo (see the comment after 
Eq.(^). This is the standard notion of decoupling of the heavy fermion. Therefore, we 
are justified in the use of Eq.(||) in our problem. 

Secondly, one could object that the rotation (^) is not a symmetry, and thus the 
subsequent manipulations are not justified. Indeed, the transformation (|1^) is not a 
symmetry of the theory, and it does change the partition function. However, the non- 
conservation of the current corresponding to the rotation ([131) shows up only at the level 
of 0(1/M^), and thus can be neglected in the limit M — » oo. The latter statement is 
simply a rephrasing of the textbook fact that a heavy quark of mass M ^ ^qcd does not 
contribute to the axial anomaly. The key difference from the case of the light quarks is 
the cancellation between the physical and regulator fields at the leading order in 

Thirdly, one could suspect that the substitution (|16D is very tricky due to the "double 
scaling" limit Kym <^ M -C Mr which is implied in such replacement. In particular, 
one may worry about a region when two densities in Eqs.(|l^,^) are close together and 
essentially form a new operator of higher dimension which has different transformation 
properties. However, this small region x ~ is under theoretical control, and the 

operator product expansion (OPE) can be used to find that the contribution coming from 
this region would be of order M"" and is thus negligible. 

Finally, we would like to discuss the problem of higher loop contributions in Eq. (^I]) . 
In the course of our transformations we have used the one-loop integration of both the 
auxiliary and Pauli-Villars fermions. Therefore, our result ( pi]) is valid up to higher loop 
contributions in correlation functions (|^) and (pT]). Moreover, in the "double scaling" limit 
AyM -C M -C Mr which is implied in our calculations, this problem is just equivalent to 
the question of higher regulator loop contributions to the conformal and axial anomalies. 
One expects that the inclusion of higher regulator loops results in the substitution of the 
one-loop anomaly —bas/{8n)G'^ by the full renormalization group invariant expression 
f3{as)/{4:as)G'^ in zero momentum correlation functions. 
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3 Witten-Veneziano formula and gluon condensate 



As an application of the relation (0), we can consider the famous Witten-Veneziano [|T2 



scheme of resolution of the U(l) problem, which is based on a large approach. We 
should note that our result (|21| ) is valid for arbitrary Nc- However, it is the large 
picture that leads to connections of physical characteristics of the 1]' with zero momentum 
correlations functions of pseudoscalar gluon currents in gluodynamics. As will be shown 
in this section, the low energy relation (0) yields a new mass relation for the t]' in terms 
of the gluon condensate in gluodynamics. 



We recall the Witten-Veneziano mass formula for the 7]' |jT2 



/,2,mJ, = -12 {-^GG , -^GG)yM , (22) 

where /^/ is the rj' residue 

(0| E (1^1^^15Q^\V') = , {^\nf-r^GG\r,') = V^f^^ml, (23) 



qi=u,d,s 



and the two-point function of the topological density in (|22D refers to the pure YM case. 
It is the latter object that is now given by Eg. ([2T|) . We thus arrive at the relation which 
is valid in the double limit 0, Nc ^ oo^ 



fv'< = Yb^fG')yM. (24) 



Our result ( P4D suggests an explanation why the correlation function (|22D is non-zero de- 
spite the fact that it vanishes to all orders of perturbation theory. In the Witten-Veneziano 
scheme, this fact was assigned to a ghost pole contribution to (|2^) . Our derivation a poste- 
riori shows that a non-vanishing value of correlation function (|2^ ) is due to a subtraction 
constant (gluon condensate) in the pseudoscalar channels, which is proportional to an 
analogous subtraction constant for the scalar channel, and is related to regulator contri- 
butions to the conformal anomaly. The ghost is simply a way to parametrize them. We 
note that on a qualitative level, this interpretation of the ghost was first discussed in ET . 



4 The 9 dependence in YM theory for 6 <^ti 

In this section we will show that the 9 dependence of the vacuum energy, topological 
density and partition function in YM theory can be found exactly for any and small 
9 <^ -K provided the relation ( pi]) takes place. More precisely, we will demonstrate that 
the Taylor expansions in the vacuum angle 9 can be re-summed exactly for these objects. 
Thus, formally, our results hold for any 9. However, it is believed that the physics is not 
analytic for 9 tc. This non-analyticity is related to the existence of distinct disconnected 
vacua which cross in energy at ~ tt. The Taylor expansion in 9 refers to a state of lowest 
energy (at small 9) out of this set, and thus can not probe other states. As a result, it can 

^ Approximate mass relations for the 77', similar to (|^, have been repeatedly discussed in the literature 
starting from |6| and . 
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be trusted only for 9 n. As will be discussed below, the existence of additional vacua is 
crucial to establish the correct periodicity in 6 with period 27r. A corresponding analysis 
TB| reveals that indeed Eqs.(P7D,(P^ (see below) correctly describe the 9 dependence for 



small 9 <ti 71. 

Let us first address the 9 dependence of the vacuum energy. The latter is defined by 
the 6'-vacuum expectation value of the trace of the momentum-energy tensor : 

= -^{^G^)e = -^({0\^G^\0) 
62 TV 62 \ TV 

+ ^l^^^dxdy{0\T{^G\0)^GG{x)^GGiy)]\0) + ..). (25) 
21 J I TT Stt Stt J / 



The three-point correlation function is now fixed by the low energy relations (^) and (pll) . 
Indeed, differentiating Eq. ([T7|) in respect to 9, using Eq.(^) and setting 9 to zero in a 
final expression, we obtain 



J dxdy {0\t[^G\0) ^GG{x) ^GG(y)} |0) = - (|)' (0|^G^|0) . (26) 

As before, the Wick T-product is implied in Eq.(^). We thus see that the three-point 
function (^) is expressed in terms of the same nonperturbative gluon condensate which 
enters Eq.(P) and (PI). Repeating the same procedure to all orders in 9 in expansion (^Sj) , 
we arrive at the following result for 9 it : 

The 9 dependence of the topological density condensate can be found analogously : 

(9|^GG|9) = -(0|i|G^|0).in(A«). (28) 

Similarly, we can calculate the 9 dependence of the partition function Zg of YM theory. 
Since log Zg generates connected n-point functions of the topological density, resumming 
the expansion of \ogZg, we obtain 

log Ze = const + z^V{0\^G^\0) cos (^9) , (29) 
32 TT \6b ) 



where V is the four- volume. A constant in Eq.(B9) can be found from the relation |22 



logZo = logZpr + ^4^(0|-^'|0) ' (30) 

62 TT 

where ZpT stands for the perturbatively defined partition function. We note that a 
perturbative contribution to the gluon condensate is absorbed in the definition of log ZpT- 
We therefore obtain 

Zg = ZpT{gl) exp |zly(0|^G2|0) cos (^^^)| . (31) 



One can see that 9 enters relations (P?]), ( PBD and ( PT] ) in the combination 6'/A'"c. Thus, 
our resuhs confirm the long standing conjecture, first arrived at within the large A^^^ 
picture, that a 9 dependence should come in the combination 9/Nc [T^. At the same 
time, these results seem to suggest a "wrong" periodicity in 9. As was mentioned above, 
this interpretation of Eqs. (|27D, p8|) and (^) would be wrong, as they are valid only 



for small values of 6* <^ vr. It can be shown ||T3[ that the correct 2tx periodicity in 9 is 
recovered when additional nondegenerate disconnected vacua, which cannot be seen in 
the naive thermodynamic limit V — >■ cxo, are taken into account. The correct periodicity 
in 9 is the property of this set of vacua as a whole, while some two of them cross in energy 



at ^ ~ TT |T3|. For small 9 <^ n these additional states do not contribute the partition 
function in the limit V — > oo, and correspondingly are not probed by the Ward identities. 
Their existence is rather revealed with an effective Lagrangian approach, see A 



corresponding analysis shows that Eqs. (27), ( pS]) and (0), as standing for small values 



of ^ <^ TT, are not in conflict with general principles of analyticity and 2ti periodicity in 9. 

Finally, we would like to comment on the scheme dependence of our results. As was 
mentioned earlier, there is a specific scheme which is implied in all our formulas. Namely, 
we use the path integral definition for all correlation functions as the derivatives with 
respect to appropriate sources (parameters I/qq or 9). This definition of zero momentum 
correlation functions implies the Wick type of the T-product symbol. It also fixes the 
rule of subtraction of the perturbative UV divergent contributions to the correlation 
functions, once such a rule is imposed for the gluon condensate. In terms of Eq.(^T|), this 
procedure means that the NSVZ low-energy theorem (^ is obtained by differentiation of 
the logarithm of Zg/ZpT with respect to I/Qq, see Eq.(|3l|). Our relation (|2l|) and its n- 
point generalizations are obtained by differentiation of the same expression \og{Zg/ Zpt) 
(^) with respect to 9. It is clear that a change of the subtraction prescription in this 
scheme will generally alter the absolute value of the condensate, but it does not change 



the general structure of these relations, see also the appendix in |]13| for more detail. 

Last, but not least, formula ( PT| ) implies a "hidden" symmetry which may be thought 
of as the origin of the new low-energy relations similar to (|2T|) . This hidden symmetry 
arises due to the appearance of the new complex parameter r = l/gl + i6'/(127r^) in 
Eq.(^T]). Indeed, Eq.(pID implies the relation (for small 9) 



{-—G^)e = const Re Miexp ( - i^9] = const Re Miexp ( -^^r ) . (32) 
Stt \ bg^ 3b J \ b J 

Thus, the nonperturbative vacuum energy depends just on this single combination of 
the coupling constant and 9 angle. This is exactly the origin of the new relation (|2T| ) 
which together with the NSVZ theorem (|9D, ([17|) can be rewritten in a holomorphic way 
||13|| . With our definition of the nonperturbative correlators as the derivatives of the 
partition function with respect to the sources r, f, the nonperturbative vacuum energy 
in Eq.(PTD is the only relevant term to differentiate. A non- holomorphic dependence on 
Qq of the perturbative part Zp^ is irrelevant in our definition for calculation of the non- 
perturbative correlation functions. At the same time, the 9 dependence appears only in 
the non-perturbative part of Zq. Therefore, the origin of our new low-energy relations can 
be understood as the interplay between the exact NSVZ theorem (|) and assumption of 
separation of the perturbative and non-perturbative contributions in Eq. (|3TD . Once this 
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assumption is made, a new complex structure emerges due to the pure non-perturbative 
origin of the 9 parameter which combines with another parameter l/g"^ into the unique 
complex combination r. 



5 Conclusions 

Low energy theorems provide a bridge between low energy and high energy physics, in- 
cluding phenomena at the boundary of the ultraviolet cut-off. The fact that the axial 
anomaly can be described in terms of regulator fermions was emphasized by Gribov 
(see also 0) long ago. In this paper we suggested using the connection between the 
axial anomaly term and the regulating Pauli-Villars fermion bilinear, given by Eq.(|I]), 
to derive new relations for zero momentum correlation functions of the topological den- 
sity. Classically, the topological density operator is a total derivative, which suggests 
that its effects can only be treated within infrared, explicitly nonperturbative methods 
such as the instanton approach. Yet, one can use another, quantum, definition of this 
object, which is given by Eq.(|I]). In contrast to the (quasi-) classical treatment, none of 
the extremely complicated problems of the latter (for example, a compactification of the 
infrared boundary, summation over different topological classes, instanton interactions, 
etc.) arise in the quantum approach. Instead, the problems are shifted to the analysis of 
the high energy behavior which is fixed by renormalizability. Therefore, we believe that a 
look at the "dark" , ultraviolet side of the topological density is very instructive as it can 
test nonperturbative infrared methods. An interesting question which could be asked in 
reference to Eq.(|l|) is whether this quantum definition of the topological density can be 
used in lattice studies of QCD and gluodynamics. 

New low energy relations, taken together with analogous formulas for the scalar chan- 
nel 0, enable us to evaluate some nonperturbative observables in gluodynamics, which 
seem to be out of the reach of other methods, within the regularization scheme based on 
the path integral. In particular, in this scheme we have found the relation between the 
topological susceptibility in YM theory, which is one of the most important characteristics 
of the YM vacuum, and the nonperturbative gluon condensate. Furthermore, within the 
same scheme we calculated the 9 dependence for arbitrary N^. and small 9 ^ it for low- 
est dimensional vacuum condensates and for the partition function in gluodynamics, and 
found that the 9 parameter always comes in the combination 9/Nc, as was expected for a 
long time. Our results can have various applications which have been already mentioned 
in the introduction. 
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